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Chiral symmetry breaking from Ginsparg- Wilson fermions 

Pilar Hernandez^^, Karl Jansen^^ and Laurent Lellouch''^ 
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We calculate the large-volume and small-mass dependences of the quark condensate in quenched QCD using 
Neuberger's operator. We find good agreement with the predictions of quenched chiral perturbation theory, 
enabling a determination of the chiral lagrangian parameter E, up to a multiplicative renormalization. 
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1. Introduction 

Spontaneous chiral symmetry breaking (S^SB) 
is fundamental to our understanding of low en- 
ergy hadronic phenomena and it is thus impor- 
tant to demonstrate quantitatively that it is a 
consequence of QCD. A natural candidate for 
such investigations is the numerical simulation 
of QCD on a spacetime lattice. SxSB, however, 
presents the lattice approach with a twofold chal- 
lenge. 

The first is that spontaneous symmetry break- 
ing does not occur in a finite volume. In QCD, a 
possible signal of SxSB is the presence of a non- 
vanishing quark condensate defined as: 



lim lim 



mim.y 



(1) 



where {qq)m.v is the condensate for finite volume 
V and mass to. The double limit in Eq. (|^) is 
rather challenging numerically! To get around 
this problem, we resort to a finite-size scaling 
analysis. This involves studying the scaling of the 
condensate with V and m as the limit of restora- 
tion of xS is approached (m ^ 0, V finite). 

Such a study requires good control over the chi- 
ral properties of the theory, which is the second 
challenge. Indeed, at finite lattice spacing, "rea- 
sonable" discretizations of fermions either break 
continuum xS explicitly or lead to extraneous 
fermion species |^]. To minimize this problem, 
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we resort to recently rediscovered [^,|| Ginsparg- 
Wilson (GW) fermions j|] which break continuum 
xS in a very mild and controlled fashion and ac- 
tually have a slightly generalized xS even at finite 
lattice spacing 

2. Light quarks on a torus 

In a large periodic box of volume V = L"^ such 
that Ft^L 3> 1, for small quark masses and assum- 
ing the standard pattern of SxSB with Nf > 2, 
the QCD partition function is dominated by the 
nearly massless pions; the system can be de- 
scribed with the first few terms of a chiral la- 
grangian 0. If, in addition, to ^ so that 

AUL ~ ^^"^ L < 1, the global mode of the chi- 
ral lagrangian field U E SU{Nf) dominates the 
partition function, leading to a regime of restora- 
tion of xS §. 

In the quenched approximation to which we re- 
strict here, topological zero modes of the Dirac 
operator induce 1/to singularities in {qq)m,v as 
TO — > 0. To subtract these contributions, we work 
in sectors of fixed topological charge. Generaliz- 
ing the line of argument given above, the parti- 
tion function Z^, in a sector of topological charge 
V, was recently evaluated Q for the quenched 
cas^. The quark condensate in sector v, pro- 
portional to the derivative of In (Z^) w.r.t. to, is 
then —Si/ = {qq)m.v.i'j such that 

% = z [I,{z)K,{z) -f I,+^{z)K,_i{z)\ + (2) 
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^We assume here for simplicity that the Nf flavors all have 
mass m. 

^The original unquenched treatment is given in 
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where z = mT,V and Iu{z), K^{z) are the modi- 
fied Bessel functions. As advertised, there is a di- 
vergence 1/m in sectors with topology. These 
terms, however, are independent of S. 

Eq. (^) summarizes the scahng of the quark 
condensate with the volume and quark mass in 
the global mode regime, as a function of only one 
non-perturbative parameter: S. Thus, by fitting 
the dependence of the finite- volume condensate in 
quark mass and volume to Monte Carlo data, we 
can extract S in a perfectly controlled manner. 

3. Ginsparg- Wilson fermions 

To perform the finite-volume scaling analysis 
outlined above, we need to be able to reach the 
chiral restoration regime without excessive fine- 
tuning and we need an index theorem to control 
the contribution of topological zero modes. Both 
these requirements are satisfied by GW fermions 
pT[ |. In particular, the leading cubic UV diver- 
gence of the condensate is known analytically for 
GW fermions and can thus be subtracted exactly. 
The resulting subtracted condensate, S^"^, how- 
ever, is still divergent: 



^suh 



(3) 



where a is the lattice spacing. The coefficients of 
the divergences are not known a priori and have 
to be determined, preferably non-perturbatively. 
For the values of m and a considered below, how- 
ever, only the quadratic divergence is important 
numerically, weaker divergences being suppressed 
by higher powers of m. A final multiplicative 
renormalization is still required to eliminate a 
residual logarithmic UV divergence in H^. 

In the present work, we use Neuberger's imple- 
mentation of GW fermions encoded in the Dirac 
operator 111 



aDN = {l + s) 



(4) 



with A ^ 1 + s — aDw where Dw is the standard 
Wilson-Dirac operator. The parameter s must 
satisfy |s| < 1. 



4. Numerical results 

We work in the quenched approximation on hy- 
percubic lattices with periodic boundary condi- 
tions for gauge and fermion fields. We choose 
/? — 5.85, which corresponds to a"-"^ ~ 1.5 GeV 
and use standard methods to obtain decor- 
related gauge-field configurations. 

To evaluate \/V^A in Eq. (^), we use a 
Chebyshev approximation, where Pn,t 

is a polynomial of degree n, which gives an expo- 
nentially converging approximation to ^1 \fx for 
X e [e, 1] jl^. The cost of a multiplication by 
-Djv is linear in n and becomes rapidly high. To 
reduce n substantially, we perform the improve- 
ments described in [y. We take s = 0.6, a value 
at which Neuberger's operator is nearly optimally 
local for /3 = 5.85 

To determine whether a gauge configuration 
belongs to the = or ±1 sectors, we compute 
the few lowest eigenvalues of D^^^D^ by minimiz- 
ing the relevant Ritz functional As pointed 
out in |l5|, it is advantageous for this computa- 
tion, as well as for the inversion of {D^j^+m){DN+ 
to), to stay in a given chiral subspace. Having 
determined the topological charge of a configura- 
tion, we then obtain the condensate of Eq. ^ in 
three volumes (8**, 10"' and 12*) by computing 



j.sub = _ ( xr' 
" V 



D 



h.c. — 



N 



l + s 



.(5) 



where the trace is taken in the chiral sector op- 
posite to that with the zero modes ||l^ and the 
gauge average is performed in a sector of fixed 
topology v. With this definition, terms ^ 1/to in 
Eq. (||) are absent]^. Three gaussian sources and 
a multimass solver jlTl were used to compute the 
trace in Eq. for seven values of m. 

We show in Fig. |^ our results for a^J^l^j.^/am 
as a function of bare quark mass. We have 15, 10 
and 7 gauge configurations on our 8'', 10"* and 12"^ 
lattices, respectively^. The solid lines are a fit of 
the data to Eqs. (||) and (||) for all volumes and 

^Though not shown explicitly in Eq. we correctly ac- 
count for the real eigenvalues of Djv at 2(1 + s)/a. 
*For the larger volumes, = configurations are rare 
while the calculation of ^"^q requires large statistics m. 
1^1 > 1 configurations, on the other hand, are rare in the 
smaller volumes. 
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masses. This fit has only two parameters, namely 
E and the coefficient of the quadratic divergence. 
We find a^E = 0.0032(4) and C2 = -0.914(8). 



Note: Related work with Neuberger's operator can 
be found in /[I^,[7^,[i^. 
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Figure 1. Mass dependence of the condensate for the 
8^ (circles), 10^ (squares) and 12* (triangles) lattices. 
The curves result from a fit to Eqs. ^ and 

Clearly, the formulae derived in quenched xPT 
give a very good description of the numerical 
data. The value of E that we extract is, in phys- 
ical units, E(/i - 1.5 GeV) = (22ll^ MeV)^, 
up to a multiplicative renormalization constant, 
which has not been computed yet for Neuberger's 
operator. The quoted error on E is purely statisti- 
cal and the statistics are rather small. Quenching 
and discretization errors, for instance, as well as 
possible contributions from higher orders in xPT 
are not accounted for. Nevertheless, the value ob- 
tained and the agreement with qxPT support the 
standard scenario of SxSB. 

We further consider the me an value of the low- 
est non-zero eigenvalue of d]^Dn in different 
topological sectors. In Random Matrix Theory 
the distributions of these eigenvalues are given 
solely in terms of E [|l8| . Our determination of E 
therefore yields predictions for the mean values. 
These can then be compared to the average val- 
ues obtained in simulation. With our 8^ results, 
we find agreement within roughly one standard 
deviation for = 1 (29 configurations) and two 
for = (41 configurations) |l|. 



9 



P. Hernandez, K. Jansen and L. Lellouch, 
hep-lat/9907022. 

H.B. Nielsen and M. Ninomiya, Phys. Lett. 
105B (1981) 219. 

P. Hasenfratz, Nucl. Phys. PS63 (1998) 53. 
H. Neuberger, Phys. Lett. B427 (1998) 353 
P.H. Ginsparg and K.G. Wilson, Phys. Rev. 
D25 (1982) 2649. 

M. Liischer, Phys. Lett. B428 (1998) 342. 
J. Gasser and H. Leutwyler, Nucl. Phys. 
B307 (1988) 763. 

J. Gasser and H. Leutwyler, Phys. Lett. 188B 
(1987) 477. 

J.C. Osborn et al, Nucl. Phys. B540 (1999) 
317; P.H. Damgaard et al, Nucl. Phys. B547 
(1999) 305. 

H. Leutwyler and A. Smilga, Phys. Rev. D46 
(1992) 5607. 

P. Hasenfratz, Nucl. Phys. B525 (1998) 401; 
P. Hasenfratz, V. Laliena and F. Nieder- 
mayer, Phys. Lett. B427 (1998) 125. 

12. H. Neuberger, Phys. Lett. B417 (1998) 141. 

13. P. Hernandez, K. Jansen and M. Liischer, 
Nucl. Phys. B552 (1999) 363. 

14. R.G. Edwards, U.M. Heller and T.R. Klassen, 
Nucl. Phys. B517 (1998) 377. 
R.G. Edwards et al, Phys. Rev. D59 (1999) 
094510. 

B. Bunk et al, DESY report (Sept. 1994); 
T. Kalkreuter and H. Simma, Comput. Phys. 
Commun. 93 (1996) 33. 

17. B. Jegerlehner, Nucl. Phys. PS63 (1998) 958. 

18. P.J. Forrester, Nucl. Phys. B402 (1993) 709; 
S.M. Nishigaki et al, Phys. Rev. D58 (1998) 
087704. 

R.G. Edwards et al, Phys. Rev. Lett. 82 
(1999) 4188. 

P.H. Damgaard, R.G. Edwards, U.M. Heller 
and R. Narayanan, hep-lat/9907016. 



10 



11 



15 



16 



19 



20 



